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Abstract 

We formulate variation principle for force-free magnetosphere of an inclined pulsar: £ + CI ■ M 
( £ , M are electromagnetic energy and angular momentum, fi is the angular velocity of a star) 
is stationary under isotopological variations of magnetic field and arbitrary variations of electric 
field. The variation principle gives the reason for existence and proves local stability of current 
singular layers along magnetic separatrices. Magnetic field lines of inclined pulsar magnetosphere 
lie on magnetic surfaces, and do have magnetic separatrices. 

In the framework of the isotopological variation principle, inclined magnetospheres are expected 
to be simple deformations of the axisymmetric pulsar magnetosphere. A singular line should 
exist on the light cylinder, where inner separatrix terminates and outer separatrix emanates. The 
electromagnetic field should have an inverse square root singularity near the singular line inside 
the inner magnetic separatrix. 

Large distance asymptotic solution is calculated, and used to estimate the pulsar power, L 
c _3 /i 2 4 for spin-dipole inclinations % 30° . 
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I. INTRODUCTION 



Pulsars are very simple and clean astrophysical objects. For all we know, pulsars are 
just magnetized spinning conductors. The only parameters of a pulsar are the distribution 
of entering magnetic field over the surface of a star and the angular velocity. For some 
observables, the magnetic dipole might be the only relevant input parameter characterizing 
the magnetic field. For these observables, pulsars form just a three-parameter family of 
objects - all pulsars with the same spin, dipole, and spin-dipole angle are equal. Yet the 
beautifully rich pulsar phenomenology remains unexplained quantitatively: given spin and 
dipole, we cannot predict radiation. 

For young pulsars, like Crab, the spin-down power is much greater then the power of 
observed electromagnetic radiation (eg. Fig. 9 in This might indicate that the magne- 
tosphere is close to a force-free state, and the magnetosphere can be calculated independently 
of radiation. Once the force-free magnetosphere is understood, one will be on firmer ground 
for predicting radiation (calculating small dissipative corrections to a pure force-free state). 
In any case, quantitative predictions of radiation require the knowledge of magnetosphere. 

The quantitative theory of force-free pulsar magnetosphere has recently appeared. Based 
on classical early insights 0, Q| , the shape, singularity structure, spin-down power, and 
MHD-realizability of an axisvmmetric pulsar (meaning pulsar with small spin-dipole inclina- 
tions) has been established [^flB]- I n particular, the spin-down power of an axisymmetric 
pulsar turned out to be 1.5 times greater than perpendicular dipole power: L « c~ 3 fi 2 Q A . 

Force-free magnetosphere of an inclined pulsar should be a doable problem too. Here we 
show that inclined pulsar magnetosphere should possess stable current /charge layer singu- 
larities. This is done using isotopological variation principle (§2 ). We also calculate the 
large-distance asymptotic of the magnetosphere, and use this asymptotic solution to give 
a crude estimate of pulsar power (§3 ). In Appendix, we discuss a mathematical problem 
of minimizing functionals by isotopological transformations. We show that isotopological 
relaxation turns separatrices into singular current layers. 
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II. ISOTOPOLOGICAL VARIATION PRINCIPLE AND SINGULARITIES 



We assume that magnetosphere is force free: 

pE + jxB = 0. (1) 

Using Maxwell equations 

<9 4 B = —V x E, p = V ■ E, j = V x B - d t E, (2) 

assuming that electromagnetic pattern rotates with velocity V = fl x r, just like it does for 
a dipole in vacuum: 

<9 t B = V x (V x B), <9fE = E • VV — V ■ VE, (3) 

and using the condition of zero tangential electric field on the surface of the star in the star 
frame, we get E = —V x B and the pulsar magnetosphere equation 

B x V x (B + V x (V x B)) =0. (4) 

The pulsar magnetosphere equation (j3J) follows from a variation principle: SS = 0, where 
S is proportional to electromagnetic action 

S = J d 3 r(B 2 - (V x B) 2 ), (5) 

and the variation of action should be calculated for isotopological perturbation of magnetic 
field 

5B = V x (<?f x B). (6) 

In other words, action should be extremized by continuously displacing magnetic field lines. 

The isotopological variation requirement is more than just a mathematical device giving 
stationary solutions. The actual time evolution of magnetic field in force-free electrodynam- 
ics is isotopological, because c^B = Vx(vxB) ([y],v = Ex B/5 2 ). Then locally, the 
variation principle requires minimizing proper magnetic energy (magnetic energy in a frame 
where the electric field is absent) by displacing the magnetic field lines. This immediately 
shows that current layers are locally stable. 

An equivalent formulation of the variation principle (jSJ |BJ) states that for stationary 
magnetosphere, the first variation of the following integral of motion of FFE (force-free 
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electrodynamics, Appendix in 



I 



d 3 r( 



B 2 + E 2 



+ ft ■ r x (E x B) ) 



(7) 
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vanishes. Here the electric field variation <5E is arbitrary and the magnetic field variation 
<5B is itopological, given by ©• 

In the framework of isotopological variation principle, the only difference between inclined 
and axisymmetric pulsar magnetospheres is the distribution of entering magnetic field lines 
over the surface of the star. It is then plausible that stationary inclined magnetospheres do 
exist. 

Another way to see that stationary inclined magnetospheres may exist, is to write the 
magnetosphere equation (jlj) as 



where A is constant along field lines: B • VA = 0. One needs to specify the scalar A on 
the surface of the star. Just like for an axisymmetric dipole [a], the right choice of A on 
the surface of the star may allow to smoothly cross the light cylinder (light cylinder is the 
surface V = 1). 

Magnetosphere equation (JSJ) shows that stationary magnetospheres possess magnetic sur- 
faces. This is because B • VA = 0, meaning that magnetic field lines lie on the surfaces 
A =const. Magnetic separatrix surfaces will then exist, between open and closed field lines 
within the light cylinder, and between oppositely directed field lines outside the light cylin- 
der. As we show in Appendix, isotopological extrema are expected to have stable singular 
current layers along magnetic separatrices. 

Moreover, in view of the variation principle, inclined magnetospheres should be just simple 
deformations of the axisymmetric magnetosphere, allowing to propose the following structure 
of the pulsar magnetosphere, based on the singularity structure found in the axisymmetric 
case: 

1. There is a singular line on the light cylinder. 

2. Inner separatrix terminates on the singular line. 

3. Outer magnetic separatrix emanates from the singular line. 




(8) 



4 



4. Separatrices are singular current /charge layers. 

5. The inner separatrix has a wedge at the singular line. 

6. In the lab frame, the electromagnetic field has an inverse square root singularity near 
the singular line inside the inner separatrix. 

This picture should be checked by an accurate numerical simulation of the inclined mag- 
netosphere. In the next section, we will just show that an asymptotic solution of the magne- 
tosphere equation exists. We use this asymptotic solution to give a rough estimate of pulsar 
power. 

III. PULSAR ASYMPTOTIC AND POWER 

At large distances from the light cylinder (for V ^> 1), the magnetosphere equation (jSJ), 
written in spherical coordinates (r, 9, 0), with Q along 9 = 0, in units ft — c — 1, is 



d e (sm 9B^) + r 2 sin 2 9d^B e = r sin 9XB r , (9) 

d r (rB^) + r 2 sin 98^ = -r\B e , (10) 

sin 2 9d r (r 3 B e ) - r 2 «9 e (sin 2 9B r ) = -r\B h (11) 

sin^^r 2 ^) + rdo(sm6B ) + rd^B^ = 0. (12) 
These equations have the following asymptotic solution: 

B T = fr- 2 + o(r- 2 ), B e = gr- 3 + o(r- 3 ), B^ = - sin 9fr~ l + o(r~ l ), (13) 



where /, g are functiuons of 9 and 0, with <fi = + r. As it should be, asymptotically one 
has E = B, and Poynting flux S oc r~ 2 . 

One can use the asymptotic solution for a crude, but "physically motivated" estimate of 
the pulsar power. Using the asymptotic solution, we calculate the pulsar power (Poynting 
flux through a distant sphere) as 

L oc J dVLf 2 sin 2 9 oc J dN(dN/dtt) sin 2 9, (14) 

where Q is the solid angle, and iV is the number of field lines. While we obviously cannot 
calculate this integral without a full solution, we can offer the following speculation. 
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The number of field lines going to infinity is equal to the number of field lines crossing 
the light cylinder. We further arbitrarily assume that (i) the number of lines crossing the 
cylinder can be approximated by the number of pure dipole (non-rotating) lines crossing the 
same cylinder, and (ii) the lines stay within the same solid angle after crossing the cylinder. 
Then (|14j) can be approximated by 

L oc J dzd(f>B 2 ± sin 2 9, (15) 

where the integral is over the light cylinder, and B± is the component of the pure-dipole 
field perpendicular to the light cylinder. Calculating the integral for aligned and orthogonal 
dipoles, arbitrarily interpolating, and normalizing to the known aligned dipole power, we 
get L »s c -3 /i 2 fi 4 (l + 0.8 sin 2 x), where x is the spin-dipole angle. Obviously, even the sign of 
the second term is questionable. This estimate, however, may serve to show that the pulsar 
power is 

L « c~V^ 4 , (16) 

for inclinations x ~ 30°. 
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APPENDIX A: ISOTOPOLOGICAL RELAXATION 

Consider the problem of minimizing magnetic energy W 

W = J d 3 rB 2 , (Al) 
by isotopological variations of magnetic field 

5B = V x (Six B). (A2) 
Performing the variation, we get the equation for the minimum energy configuration, 

B x V x B = 0. (A3) 
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FIG. 1: Initial ip isolines, and ip isolines some time after the incompressible relaxation (|A6|) was 
turned on. 

This model is a simplified version of the full pulsar case. 

Since W is positive, any initial configuration can be relaxed into a minimal energy con- 
figuration by a series of isotopological deformations. We claim that these deformations turn 
separatrices into current layers. 

A simple rigorous proof can be given after some further simplifications. Consider two- 
dimensional magnetic fields B = (—d y tp,d x ijj,B), where ip and B are functions of x and y 
but do not depend on z. We also assume that minimizing transformations are z- independent. 
If B 2 3> (V^) 2 , the field first relaxes to a state of nearly uniform B by compressible defor- 
mations in the xy plane. Then, the magnetic energy of the plane field 

W = J d 2 r{V^) 2 , (A4) 

will be minimized by incompressible xy motions, that is by 

S^ = {S X M, (A5) 

where {x, = d x xd y ip — d y xd x il> is 2D Jacobian. 

One can implement an explicit minimization. For example, 

^ = {X,^}, X=W,VM (A6) 

is an incompressible motion giving W < 0. We solved (|A6|) numerically. Fig.l shows 
how an X-point splits into two Y-points. The inner part of the separatrix is obviously a 



singular current layer. The outer part is also a singular current layer. To see this, note that 
§ dl\ V^| _1 along each -0-isoline is conserved by incompressible motions, while the inner and 
outer lengths § dl of the relaxed separatrix differ. 

Fig. 2 gives a simple way to understand the singularity (and, in fact, can be developed 
into a rigorous proof). If most of energy is in the inner part, the relaxation transforms the 
separatrix into two circles. If most of energy is in the outer part, the relaxation transforms 
the separatrix into one circle. Hence, a generic initial state gives a final state with two 
Y-points. 
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FIG. 2: Initial ip isolines, and tp isolines after incompressible relaxation. 
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